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1. introduction 

The labor of computing terms or smns of terms in a binomial ex- 
pansion can be greatly reduced by the use of certain continuous func- 
tions approximating to the values desired. Although these functions 
are treated in many texts and papers, the simple transformations and 
substitutions required to obtain formulae directly applicable to com- 
putations pertaining to binomial series are in general left for the stu- 
dent to make. 

This paper brings together the most suitable of these functions and 
illustrates their use in order to meet the needs of investigators wishing 
to obtain numerical results as readily as possible. 

2. THE symmetrical BINOMIAL 

The xth. term counted from the middle one in the expansion of the 
symmetrical binomial* 

where the number of the term is 

r = l, 2, 3, 4 ... (2) 

is approximately given by the familiar exponential expression (Yule> 
1911, p. 300; Merriman and Woodward, 1900, pp. 482-483) 

* When n is even, there are an odd number (n +1) of terms, and corresponding to the middle term 
which is also the maximum, x has the value zero. 
When n is odd, there are an even number (n +1) of terms, the two central ones are equal, and the 

corresponding values of x are, respectively -{ — and . 

2 2 



In both 



r = -+H-x 
2 



n+1 1 

x=r ■ . 

2 2 
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where 
and 



x' x' 

y = ^=-e~^ = ^^e~^ (3) 



ff = the standard deviation = ^ Vn (4) 

e = 2 . 7183, the Naperian base of logarithms. 



With the aid of a tabulation of this exponential function (Davenport, 
1899, p. 21 and 54; Fowle, 1914, pp. 48-54) its value can be quickly 
found when n and x are given. The difference between this function 
and the corresponding term of the binomial series decreases as n in- 
creases, but the approximation is close even when n is as small as 8, 
as shown by Table I. 

TABLE I 

VALUES OF TERMS IN THE BINOMIAL EX- 
PANSION COMPARED WITH THE COR- 
RESPONDING EXPONENTIAL FUNCTION 



r 


X 


V 


Terms of the 
expansion 


1 


-4 


.0052 


.0039 


2 


-3 


.0313 


.0313 


3 


-2 


.1038 


.1092 


4 


-1 


.2198 


.2184 


5 


-0 


.2822 


.2730 



To compute the value of y, substitute 8 for n in the general expres- 
sion 

2X' 

(5) 



.798 -- 
y= ^ 



Vn 
obtained from equations (3) and (4), thus obtaining 



y=.2822e 4. 



(6) 



To find y for a given value of x (x is negative for terms to the left of 
the middle one), say, x=—2, substitute — 3 for x in equation (6) and 

-- 1 

find the value of e * =e-i = - from table 18 (Fowle, 1914, p. 49). 

e 

Thus 



2/ = (.2822) (.3679) = .1038. 



(7) 
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To use Davenport's (1899) table (p. 54), find from equation (4) 

ff = iV8=V2. (8) 

Thus for a; =-2 



-=--7= = -V2 = -1.414 



y 



1- =e ^' = .3753- .14(. 3753- .3246) = .368 

I/O 

and as before 

2/ = (.2822) (.368) = .1038. 



(9) 
(10) 

(11) 



If n equals the odd number, 9, substitution in equation (5) gives 

2/=.266e 9 , (13) 



from which the entries in column 3 of Table II were computed. 



TABLE II 

VALUES OF TERMS IN THE BINOMIAL EX- 
PANSION COMPARED WITH THE COR- 
RESPONDING EXPONENTIAL FUNCTION 



r 


X 


V 


Terms of the 
expansion 

8+J)' 


1 


-4.5 


0.0030 


0.0020 


2 


-3.5 


0.0175 


0.0176 


3 


-2.5 


0.0663 


0.0704 


4 


-1.5 


0.1622 


0.1640 


5 


-.5 


0.2516 


0.2460 


6 


.5 


0.2516 


0.2460 



A continuous function approximating to the terms of a series is 
especially useful in computing the sum of a given number of terms, 
since another continuous function approximating to this sum results 
from integration. Thus the sum of any number {x) of terms counted 
from the middle or zero'th term of the binomial expansion is given 
approximately by the integral 



rx . rx .± lY f^ -^ 

/ ydx = -^\ e 2>^da; = -v/— / e » dx 
Jo "^l-KaJ" ^imjo 



where y is given by equation (3), or (5). 



(14) 
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/2 
Substitution of a new variable t = x\j— in equation (14) gives the 

^n 

well-known "probability integral"* 



/ ydx = -7= / e-^'dt 

Jo V-^Jo 



(15) 



whose values are tabulated in various books on statistics and least 
squares (for example, Yule, 1911, p. 306; Davenport, 1899, p. 55; 
Broggi, 1911, pp. 354-356; Fowle, 1914, p. 56; Merriman, 1903, p. 220). 

When n is even, if x has the value ±|, the integral (14) or (15) gives 
^ of the middle term, for a; = ±(l+^) the integral gives the value of 
i the middle term, plus the next term, and so on. 

When n is odd, if x has the value ±1, the integral (14) or (15) gives 
one of the two central terms, for a; = ±2 the integral gives one of the 
two central terms plus that of an adjacent term, and so on. 

To find the sum of the first r terms [counting the 1st term as one, 
see equations (1) and (2)], find the sum from the middle term to, but 
not including the r th term, and subtract the result from |, which is the 

n+1 
sum of the first terms. That is, for odd or even values of n 



■'<-'¥) 



and the sum of first r terms is 



where 



(2 _ n-2r+l 



(16) 



(17) 



V^ ■ (18) 



Table III gives the sum of portions of several binomial expansions, 

* To determine the probability of obtaining any proportion of an equal number of each of two different 
things, put 



where is the difference between the true proportion, J and the observed proportion, then 



This agrees with the result derived by Hatai (1910) from the theorem of Bayes for determining the 
probability that the true proportion differs by 6 from the observed value. 
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taken from a biological paper (Sumner, 1910, p. 330), and the approx- 
imation to these sums given by the "probability integral." 



TABLE III 

SUMS OF TERMS IN CERTAIN BINOMIAL EXPANSIONS AND CORRESPONDING 
VALUES OBTAINED FROM THE PROBABILITY INTEGRAL 



n 


r 


«-2r+l 


Vm 


£ equals 

»-2r+l 

V2n 


¥> 


-A/> 


'[-*/>] 


Sum of 
terms of 
binomial 


10 


2 


7 


4.472 


1.566 


.9732 


.0268 


.0134 


.0193 


11 


1 


10 


4.691 


2.133 


.9974 


.0026 


.0013 


.0005* 


11 


3 


6 


4.691 


1.280 


.9297 


.0703 


.0351 


.0327 


11 


4 


4 


4.691 


.853 


.7725 


.2275 


.1137 


.1111 


18 


6 


7 


6.000 


1.167 


.9011 


.0989 


.0494 


.0588 


19 


4 


12 


6.160 


1.948 


.9941 


.0059 


.0029 


.0022 


19 


6 


8 


6.160 


1.298 


.9336 


.0664 


.0332 


.0318 


32 


7 


19 


8.000 


2.375 


.9992 


.0008 


.0004 


.00041 


56 


14 


29 


10.580 


2.745 


.9999 


.0001 


.00005 


.00005 



* The approximation is closer as r increases, and for computing one or two terms the binomial series 
itself should be used, instead of the integral. 

To illustrate the process, determine the sum of the first 7 terms in 
the expansion 

(i-r=©--©''-if ©-- • • • 

r = l, 2, 3, . . . (19) 



Here, 



and t = 



n=32, r = 7, n-2r+l = 19 



(20) 



19 



V2.32 



= 2.375 (Table III, line 8). The corresponding value 



2 n 

of the integral -7=/ e-*'dt (Merriman, 1903, p. 220), is .9992, and 

1(1 — .9992) = .0004 the approximate value, which agrees well with 
.00041, the sum of the first 7 terms of the series (equation 19). 

To use Table IV (Davenport, 1899, p. 55) the above expression for 
t must be replaced by 

n+1 
X 2 



•— r 



ff Wn 



w-2r+l 
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and the tabular entry must be subtracted from J. If n = 32, and r = 7, 
.^32-14+1^ 19^ ^3 3^^ 
ff V32 V32 

and the corresponding tabular entry is .4996. The approximation to 
the sum is therefore .5000- .4996= .0004. 

3. THE ASYMMETRICAL OR SKEW BINOMIAL 

A function bearing the same relation to the skew binomial 

r = l, 2, 3, 

as the normal probability function bears to the symmetrical binomial 
was first found by E. L. Deforest (1882-3). In a review of Deforest's 
researches, Shinkishi Hatai (1910) presents a derivation of Deforest's 
(1882-3) formula 






kV2irb\ «&/ * (22) 

Where y is the approximate value of the a^h term in the binomial 
series, the mean term r = (ng+1) being the origin, 

A; = H — + f- . 

12a^fe 2S8{aW ^ ^23) 

a= — , 0=112 

till and ixz are respectively the second and third moments about the 
mean. These moments* of the binomial series are (Pearson, 1906-7, 

2/' x' 

*By definition, tlie second moment /«j = , where / is the frequency of a deviation x, and 

N 
N =Xf'= the sum of all the frequencies. If the interval between successive terms eqiials L, instead 
of unity, as assumed in the preceeding pages, 

lit"— , 

N 

since the deviation is xL. For a binomial series, the corresponding standard deviation becomes 



If the total range, nL is denoted by R, the standard deviation becomes 



!, the standard deviation be 

"C^)s-(V?)»- 
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p. 172; Greenwood, 1913, p. 71) 

lJ.3=npq(p-_g) } 

To illustrate the use of Deforest's rule, compute the rth term of the 
expansion ( . 9 + . 1)^". By definition 

p=.9, q=.l, n=20, 20X. 1+1=3 
mean term or origin, 

M2=(20)(.9)(.l) = 1.8 

M3 = (20)(.9)(.l)(.9-.1) = 1.44 

= ^^^=2.6,6 = 1.8, o«& = ll. 26 
1.44 



jS; = H 1 H 1 =1.0076. 

(12)(2. 5)^1 .8) (288)[(2. 6)^(1 .8)P 



Therefore 



■H')' 



.10.25 



(26) 



where x = {r— mean term) = r — 3. 
To obtain the third term put x=0, then y= .295. 
Ta obtain the 6th term, put a; = 3, then 

log .295= 9.46982-10 

10 . 25 log (- ] = (10 . 25) ( . 221858) = 2 . 27404 
-7.5 log e = (-7.5)(.43429) = -3.25717 



log 2/= 8.48669-10 
y= .03067. 

As shown by Hatai (1910), Deforest's formula is equivalent to 
Pearson's Type III curve whose equation (referred to the mode as 
origin) is (Pearson 1914, p. xvl) 



y=y<fi \}+—) 



ml ' (26) 

In computing terms of a binomial expansion by means of equation 
(26), the constants have the values (Pearson, 1906-7) 

4 , 2 

«i = 1, v = 



4 p — q 

npq — 
n 
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\m+l 
But a more convenient form 

y=yoe '-(1+-) 

ai \ ai^ (.^1) 

for numerical work is given by Pearson (1914, pp. xlv and 37). These 
two forms will be equivalent if pi = m and ai = — . Substituting the 

V 

values of m and v gives 



p,- 


4 


-1 = 


_4ff2(n+l)- 
n-i<T^ 


-n 


•Tl- 


1 4 






ff« n 








«! = 


.Pi(p- 
2 


s) 







and 



where <r* = npg. In formula (27) the origin is at the mode, hence 
;i;=r— mode. 

But for Pearson's Type III curve the relation of the mode to mean 
is (Elderton, 1906, p. 65) 

Mode = mean — ^ — 

which reduces to 

Mode =nq+l — - — - for the binomial series. Therefore in formula 
2 



(27) 



a;=r— ng— 1-"- — 



2 ■ (28) 



Also the simpler expression 

1 .399 1 



V27r<r <r ^2'Knpq 

can be used for yo (Merriman and Woodward, 1900, 482). Table IV 
gives the terms of the expansion of (.9+.1)*" computed from the 
series, and from Pearson's Type III curve. 
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TABLE IV 

TERMS OF THE EXPANSION (.G+.l)"" AND THE CORRESPONDING VALUES OF 

PEARSON'S TYPE III FUNCTIONS 



r 


X 




-..^ 


Ko» 


V 


Value of 
terms, com- 
puted from 
the expansion 


1 


—1.6 


—.390 


.046 


1.001 


.1002 


.122 


2 


—0.6 


—.146 


.0052 


1.420 


.263 


.270 


3 











1.473 


.297 


.285 


4 


1.4 


.342 


.020 


1.268 


.185 


.190 


5 


2.4 


.585 


.052 


2.939 


.087 


.090 


6 


3.4 


.830 , 


.098 


2.466 


.029 


.032 


7 


4.4 


1.073 


.149 


3.944 


.009 


.009 


8 


5.4 


1.318 


.208 


3.340 


.002 


.002 


9 


6.4 


1.560 


.269 


4.713 


.0005 


.0004 


10 


7.4 


1.805 


.336 


4.025 


.0002 


.0001 



To illustrate the use of Pearson's formula (equ. 27) compute y when 
r = 5. Thenn=20, p=.9, (?=.l, 

ff* = 20X.9X.l = 1.8, ff = 1.342 
4ff2 = 7.2 



^ 7.2X21-20 131.2 ,. „. 10.25(.9-.l) ., 

Pi = = = 10.25, ai= i^ i-=4.1 

20-7.2 12.8 2 



and 



y,= ^^ = .2973 and X = -. 
1.342 oi 



To find log loy, subtract Pi times the value of 
{logio(l+X)-XLogioe} 

found in Pearsons' Table XXVI (1914, p. 37), from logw j/o. 

2 4 
For example, if r = 5, a; = 5-3 + .4=2.4, X=-^ = .585 

{logio a+X)-X Logic e} = .052 
(10.25) (.052) = .534 

log 2/0= log .2973 = 1.473 
log2/ = 2.939 
3/ =.087 
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The sum of any number (r—s) of terms measured from the mean is 

fydx = j/o re-^'"'(l + -Tdx 

Jo Jo \ Ui' 



(29) 



where x=r—s+- — -: 
2 



but tables of this integral are not available. 

Edgeworth's (1906) theory of skew variation based on an extension 
of the methods by which the norman law was derived yields much more 
convenient formulae than Pearson's more general theory. ToUey 
(1916) concluded that Edgeworth's (1906) approximate rule (see 
equation (30)) would suffice for climatic data, and has published con- 
venient tables for its application (Tolley 1916, pp. 640-1). 

The total area of the generalized curve between the abscissa 
xi and X2 is 

^'L'i7ky^'''-[w^HTh)i$-')]z: (30, 

(Tolley, 1916, p. 638), where 

<r= standard deviation, 

x = distance (abscissa) from mean, x^ 
/ij= third moment about the mean 
N = total number of observations. 

For the binomial series (equation 21) on page 611 

p-q 



ff=npq, H2=npq (p—q), k = 



^ npq 



If y=- , — , equation (30) becomes 
■V2<r 



^=r^'-'^'-fe"'^-'>5]i" (31) 



If j/2 equals oo , A = A„ area to the right of j/i. 
If 2/1 equals — oo , A =^lj area to the left of y^. 



616 American Statistical Association [84 

That is, 

^^-^S.'-'"-\i-/-^''^-»l] (33) 

but 

v;/-/-'''*-H'-v;i, '^'«]''"<'' (37) 

and the "probability integrals" in the second members are given in 
various tables, e. g., Broggi (1911, p. 354), Czuber (1914, p. 437), 
Fowle (1914, p. 56), Merriman (1903, p. 220) and Pearson (1914, pp. 
9-10). The exponential e~^' has also been tabulated by various 
authors, e. g., Fowle (1914, p. 48), or may be obtained from a table of 
the "probability integral" since 

fy+Ay . n ry+Av o fy—Ay 



fV+Ay n fy+Ay 9 fy-Au ^ 

/ e-''dt ,_{^ e-^dt--^ e-"dt) 

Jv-Ay V^/V^JO V^JO ) 



2Ay 2 V 2^y 



(38) 



For example suppose we wish any term or the sum of any number of 
terms of the series 

(.9+.l)2«=.122+.270+.285+.190+.090+.032 
'•= 1, 2, 3, 4, 5, 6, 

+ .009+.002+.0004+.0001, (39) 

7, 8, 9, 10. 

Then 

p=.9, g=.l, n=20 
mean=ng+l =2.0+1 = 3.0 



(r = V(20)(.9)(.l) =1.342 
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fc=^-^^ =.596, - = .0993 
1.342 6 

= .5268a;, and equations (33), (35), and (37) become 



V2<7 

^L = I [l + ^ /V'd<l - . 0396 e-^^i2y,^ - 1), (40) 

for x>0, and 

AL=^\l--^r'e-''dt'\ - .0396 e-^''{2y,^-l), (41) 

for x<0, where the absolute values of y are to be substituted. 
Let ''^=0, then 2/2 = and Equation (40) becomes 
Ai = H-. 0396 =.5396. 

Let a; = 1.5, then 2/2=. 7902, yi^= .Q2U, e-"''=.5404 and Equation 
(40) gives 

i:^^-. 0053 =.8628. 
2 

Let a; = — 1 . 5, then j/2 = — . 7902 and equation (41) gives 

^~-^^^^ -. 0053=. 1266. 
2 

The table of values of Al published by ToUey (1916, pp. 640-641) 
can be used to facilitate the work, as illustrated by the following com- 
putation of the values already obtained. 

As before, fc has the constant value . 596 for all the terms of the series. 



Suppose a;=0, then from the table 
Tabular entry, Al, 
.53989, 
.52660, 


k 
.600, 
.400, 




Diff. .01329 
Diff. .00027 


.200 

. 004=. 600- 


-.596 



. 53962 = Interpolated value of Al- 

Suppose a; = 1 . 5, then - = 1.117 and from the table 
<r 

X 

Tabular entry, Al value of - 



618 

when fc = . 6 


Amt 

.85976 
.86829 


Diff. 
Diff. 


.00853 
.00289 
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1.10 
1.15 



.05 

.017 = 1.117-1.10 



.86265=Ai, when fc = 6.00 and - = 1.117 

<r 

Diff. . 00003 correction for reducing to A; = . 596. 



.8627 Final interpolated value of .4 X, 



Suppose a;= — 1 .5, then - = — 1 . 117 and from the table 



Tabular entry, 

when fc= 0.600 

.13109 
.11843 


Ai 


value of - 
<r 

-1.10 
-1.15 


Diff. .01266 
Diff. .0042 


.05 

.017=-1. 10+1. 117 



.12687= Ax, when fc = .600 and - = 1.117 
.00003 correction for reducing to fc= .596 



. 1269 Final interpolated value of Ai,. 

To compare these results with the exact values given by the series, 

we must remember that x is counted from the mean term np+1 =3, 

as origin and therefore in general (See equation 21), r = 3+x. 

285 
The first two terms plus i the mean term .122+.270+ ^ = .5345 

corresponds to x=0, and agrees well with .5396 obtained from the 
formula. When a; = — 1 . 5 the sum is simply the first term, . 122. 
The formula gives . 127. When a; = 1 . 5, the sum is the first 4 terms 
.122+.270+.285+.190=.867. The formula gives .863. 

In the case of an extremely asymmetrical binomial in which -^ is 

pn 

very small, the computation of the general term can be most readily 
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made by means of Poisson's exponential expansion. (Pearson 1914, 

p. IxxAdi and pp. 113-21.) The smaller the ratio — , the more closely 

pn 

the terms of the binomial series 

1 ■ A 



approximate to those of the exponential series 



e-^'C 



e-"* 1+m-f 



m^ 



m' 



r = 



1-2 1-2.3 
1, 2, 3, 



) 



where m=nq. The expression 



'»-"(5) 



for the (a;+l)st or rth term is 



tabulated for m=.l to 15, and x = to 37 by Pearson (1914, pp. 
113-21), thus after computing the constant m for a given case, the 
successive terms of the series can be read off directly from his 
table. For example, in the expansion (.9+.!)^' already considered 

m=ng = (20)(.l)=2 and the rth term equals e~H ) whose 

\(r — 1)!/ 

values are entered opposite those of the series in Table V. 



TABLE V 
TERMS OF THE BINOMIAL SERIES (.9H-.1)», 
AND THE CORRESPONDING VALUES OF 
POISSON'S EXPONENTIAL EXPANSION 



r 


r— 1 equals x 


e-^/2r-l\ 


Terms of 
(.9+.l)» 


1 





.135 


.122 


2 


1 


.271 


.270 


3 


2 


.271 


.285 


4 


3 


.180 


.190 


5 


4 


.090 


.090 


e 


5 


.036 


.032 


7 


6 


.012 


.009 


8 


7 


.003 


.002 


9 


8 


.0008 


.0004 


10 


9 


.0002 


.0001 
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